Finite— size scaling of the helicity modulus of the two— dimensional 0(3) model 
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Using Monte Carlo methods, we compute the finite- 
size scaling function of the helicity modulus T of the two- 
dimensional 0(3) model and compare it to the low tempera- 
ture expansion prediction. From this, we estimate the range 
of validity for the leading terms of the low temperature expan- 
sion of the finite-size scaling function and for the low temper- 
ature expansion of the correlation length. Our results strongly 
suggest that a Kosterlitz-Thouless transition at a tempera- 
ture T > is extremely unlikely in this model. 



We report a computation of the finite-size scaling func- 
tion of the helicity modulus of the two-dimensional 0(3) 
model on the lattice in the presence of periodic boundary 
conditions. The lattice hamiltonian of the 0(3) model is 
defined as follows: 



H 



(1) 



where the sum runs over nearest neighbors, the pseu- 
dospin s = (s x , s y , s z ) with s\ + s 2 y + s 2 z = 1, and J > 
sets the energy scale. The continuum model and its lat- 
tice version (||) are supposed to describe the same long 
wave length physics, i.e. they exhibit asymptotic freei 
or in other words, their critical temperature is zeroEJ 
This is in accord with the Mermin-Wagner theoren 
which forbids a transition to a state of long range order in 
a two-dimensional system of O(N) symmetry with N> 1 
at a temperature T > 0. However, a Kosterlitz-Thouless 
like transition is still allowed as is well known for the 
two-dimensional 0(2) modeE. Indeed, there have been 
arguments in favor of a Kosterlitzi-Thoulcss transition 
in the two-dimensional 0(3) modelEI which have recently 
been subject to various Monte Carlo investigationsQtl fo- 
cusing on the low temperature behavior of the correla- 
tion length. An easily accessible quantity signalling the 
Kosterlitz-Thouless transition is the helicity modulus. If 
a Kosterlitz-Thouless transition occurs in the model, on 
an infinite lattice the helicity modulus exhibits a jump at 
the transition temperature T c > from zero to a fini 
value and remains nonzero for temperatures T < T c 
However, our result for the finite-size scaling function 
suggests that the helicity modulus vanishes at all tem- 
peratures T > on infinite lattices, thus a Kosterlitz- 
Thouless transition is very unlikely to occur in the two- 
dimensional 0(3) model. 

We compare our result for the finite-size scaling func- 
tion to the low temperature expansion prediction of 



Brezin et al.0 obtained for the 0(3) nonlinear a model 
and given by 



T(r,L) = \_ ( £t(T) 
T 2tt V n L 



Inln^T] -Oil i. (2) 



Here £r de 
behaves aa 



jotes a correlation length of the system which 



frCT) - C T (27r/r)- 1 exp(27r/T). 



(3) 



The prefactor Cx cannot be determined within perturba- 
tion theory. Note that the scaling form (||) is only valid 
in the limit £ r (T) > L. 

Numerically we determine the scaling function as fol- 
lows. We compute the helicity modulus at various 
temperatures and for different lattice sizes and plot 
T(T,L)/T versus \n(LC r / £ r (T)) . In the range where 
the low temperature expansion predictions Eq.(||) and 
Eq. (||) hold, we expect our scaled data to collapse onto a 
single curve. Thus, we are able to estimate the range in 
the variable z — \n(LCr / £,r (T)) where the leading be- 
havior (||) sets in and the temperature range where the 
low temperature expression (|J) is valid. For the latter let 
us assume that the low temperature expansion (Q) were 
not valid above a certain temperature T + say. Then we 
would see deviations from the scaling curve for values of 
z larger than \n(LC r / £ r (T + )) . 

For our Monte Carlo simulations, we use the hamil- 
tonian (Q) to compute the helicity modulus of the 0(3) 
model onLxL lattices with L = 1CL40, 60, 100, 200 em- 
ploying Wolff's 1-cluster algorithmO. In order to avoid 
boundary effects we apply periodic boundary conditions. 
Though the leading terms of the scaling expression (^) 
were derived for fixed spin boundary conditions, the scal- 
ing form (^|) remains valid for periodic boundary condi- 
tions as well. The definition of the helicity modulus along 
the space direction /i in the presence of periodic boundary 
conditions is (for the derivation follow the steps outlined 
in RefsO) 



»t+A< 



2T ( 



(4) 



[l,m] 



I m m I 



where s l and s m denote two different components of the 
spin s. The summation is over all possible pairs of spin 



1 



components, denoted by J2u,m]i anc ^ over an lattice sites, 
denoted by Ylv The symbol i + /j, means the adjacent 
lattice site of i in the space direction [i. In the fol- 
lowing we will write T instead of T M as our lattice is 



isotropic. Our definition of the helicity modulus 
sures limT~*o Y(T, L) = 1 which agrees with Eq.(b) 
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FIG. 1. The helicity modulus T(T,L) as a function of 
temperature T for various lattices L x L. The error bars are 
smaller than the symbols. 

Fig.[j] shows our data of the helicity modulus as a func- 
tion of temperature for various lattice sizes. From this 
data it is impossible to draw conclusions about the behav- 
ior of the helicity modulus of the 0(3) model on an infi- 
nite lattice. For this reason we have examined the finite 
size scaling function in the manner described above. 

In Fig.| we plot T(T, L)/T versus z = InL - 2n/T + 
ln(27r/T) using the data of Fig.|l|. The data collapse onto 
a single curve with very little scatter for values z < — 5, 
only for z > -5 do the data for T(T, 10)/T begin to 
deviate from the single scaling curve (cf. the inset of 
Fig.^j), signalling the break down of the expression (^). 
The corresponding temperature is T ~ 0.657, i.e. for 
temperatures T smaller than 0.657 the low temperature 
expansion result (^) should be a good description of the 
temperature dependence of the correlation length. This 
temneputufie is larger than the value 0.37, obtained in 
RefsE'EJiL3 where the correlation length was computed 
directly via its second moment definition, and obtained 
by different methods in RefsO. 

So far we have only shown that our data of the helic- 
ity modulus computed at various temperatures and for 
different lattice sizes obey scaling. Now we wish to com- 
pare the scaled data in Fig.|| with the low temperature 
prediction Eq.(0). To this end we fit the expression 



f(z) = ai(a 2 - z + ln(a 2 - z)), 



(5) 
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FIG. 2. Scaling plot T(T,L)/T versus z = InL — 2n/T 
+ ln(27r/T). The error bars are smaller than the symbols. 



to our scaled data treating a\ and a-i = In Cx as fit pa- 
rameters. In principle we should add a constant 03 to 
f{z) to allow for the 0(1) term in Eq.(||), however 03 
turned out to be zero within error bars, therefore we set 
03 = for our fits. 
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TABLE I 
data shown in Fig|| 
quality of the fit 



Fitting results of expression (^) to the scaled 
X 2 is per degree of freedom and Q is the 



Table I contains our fitting results. We varied the num- 
ber of data points in the fit to check the stability of our 
fitting results and to estimate the range of validity of 
the low temperature result (Q) in terms of the variable z. 
We estimate this raage from the quality of fit Q (see e.g. 
Eq.(6.2.18) of Refia); the closer Q is to 1, the more reli- 
able are the values obtained for the fit parameters a\ and 
a-i- The value of the parameter a\ agrees rather well with 
the predicted pKalue of l/(27r) = 0.15915..., compare the 
results of Ref.li2 where the best value ai = 0.162(4) was 
obtained. Fig.|| shows the fitting results of the fit of the 
form (H) to our scaled data in the range —56 < z < —35, 
i.e. including 70 data points. 

Since the quality of the fit decreases by an order of 
magnitude when about 100 data points are included in 
the fit (cf. Table ||), we are on the safe side assuming 
that formula (0) is valid for values of z in the interval 
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fit: f(z)=a,(a 2 -z+ln(a 2 -z)) 
a,=0.15912(5) 
a 2 =-2.282(16) 
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FIG. 3. Fit of the expression (g) to the scaled data of 
Fig.^ in the range —56 < z < —35 (70 data points). The 
error bars are smaller than the symbols. 



( oo, 35] (cf. Fig.||), i.e. expression 
peratures and lattice sizes that fulfill 



In 



T 2ir\ 

exp —| > 35. 



2nL 



T J 



holds for tern- 



(0) 



Our results suggest that the correlation length entering 
the finite-size scaling form (j^) depends on the temper- 
ature according to Eq.(]3) with C'r — 0.102. The pref- 
actor G'fp-iis about eight times larger than the prefactor 
0.0125Lrll3 obtained from,th,e second moment definition of 
the correlation length £0'E3. Suck a comparatively large 
prefactor was also found in Rcf.EJ. The difference in the 
prefactors is no surprise to us since we can consider the 
scaling expression (^J) as providing an alternative defini- 
tion of the correlation length with a prefactor different 
from that of the second moment definition. 



differ somewhat from the predicted value of 1/(2tv). Note 
that as far as the quality of the fit Q (cf. Table ||) is 
concerned, Eq.(^) is as good a description of the scaled 
data in the range -56 < z < -35 as Eq(|). However, 
the fit value a\ «_C|.1626 does not agree as well with 
the predicted valudiJ 1 / (2tt) . The inclusion of the In In 
correction term yields much better agreement between 
the value of a\ and the universal value l/(27r) predicted 
by the low temperature expansion (Q). 

Let us conclude with the discussion of the possibility of 
a Kosterlitz-Thouless transition in the two-dimensional 
0(3) model. Such a,transition has been suggested by 
Seiler and PatrasciouD and its existence, has been subject 
to various Monte Carlo investigationsQB. Since we ob- 
serve finite-size scaling for the helicity modulus we con- 
clude from our scaled data (cf. FigM that the helicity 
modulus in the limit of an infinite lattice vanishes at all 
temperatures. Namely, by fixing the temperature and in- 
creasing the lattice size we move along the scaling curve 
in the direction of increasing z until the helicity modulus 
vanishes. The existence of a Kosterlitz-Thouless transi- 
tion, however, would imply a finite value for the helicity 
modulus below a critical temperature even for an infinite 
latticeO. Our confirmation of the validity of the low tem- 
perature expansion together with the strong indication of 
our scaled data that the helicity modulus on an infinite 
lattice vanishes, make a Kosterlitz-Thouless transition in 
the two-dimensional 0(3) model extremely unlikely. 

In conclusion we have computed numerically the finite- 
size scaling function of the helicity modulus of the two- 
dimensional 0(3) model and have found very good agree- 
ment with the results of the low temperature expansion. 
We have estimated the ranges of validity of the low tem- 
perature expansion in the corresponding variables. Our 
results suggest that there is no Kosterlitz-Thouless tran- 
sition in this model. 

N.S. would like to thank S. Burnett for valuable dis- 



cussions. 
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TABLE II. Fitting results of the expression ([?]) to the 
scaled data shown in Fig.^. \ 2 i s P er degree of freedom and 
Q is the quality of the fit. 



We have also fitted the expression 
f(z) = a 1 {a 2 - z), 



(7) 



to our scaled data (the constant a 3 has been absorbed 
into 02), i.e. neglecting the In In term in Eq.(||). The 
results are given in Table |J. We find values for aq that 
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